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Lattice studies of charmonium systems have indicated that in a deconfined glu- 
onic plasma ground state charmonia survive as bound states upto temperatures 
~ 2T C . After surveying the methodologies used in reaching these results, we exam- 
ine the behavior of these systems when the bound state is in motion with respect 
to the heatbath frame. We find that the finite momenta charmonia show medium 
modifications when the medium is deconfined; in particular, a modification of the 
energy-momentum dispersion relation is indicated. 



The study of strongly interacting matter at high temperatures is impor- 
tant both for understanding QCD and for cosmological applications, and is 
a very active field with experimental data available from relativistic heavy 
ion collisions. Charmonium systems have been of great interest to the rela- 
tivistic heavy ion community ever since Matsui and Satz, based on potential 
model studies, predicted that charmonium bound states cannot survive in a 
deconfined plasma, and suggested J/ijj as a probe of the confinement status 
of the plasma 1 . Later, more detailed potential model studies suggested that 
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all charmonium bound states dissolve before or pretty soon after deconfine- 
ment 2 . However, a direct lattice investigation of the hadronic correlators 
have revealed that the ground state charmonia (r] c and J/tp) survive the 
deconfinement transition in pure gauge theory 3,4,5 , and may dissolve only 
at temperatures ~ 2T C . 

Both for understanding the mechanism of bound state dissolution in 
the plasma and for study of charmonium abundance in experiments, it is 
important to also know the fate of these bound states when they not at rest 
but moving with respect to the heatbath frame. The lack of covariance in 
finite temperature system means the pole of the inverse propagator, 



can be a generic function of spatial momentum p, leading possibly to a 
modification of the energy-momentum dispersion relation and a momen- 
tum dependent decay width. Such modifications have been discussed in 
different models for pions in hadron gas 6 , and for decay width of J/tjj 
in hadron gas 7 . Here, we present a preliminary report of a direct lattice 
investigation into the fate of the ground state charmonia with a finite spa- 
tial momentum (in the heatbath frame). In this report, we will restrict 
ourselves to comparatively moderate temperatures < 1.5T C . 

Direct lattice investigations of the fate of the charmonium states involve 
analyzing the Matsubara correlators, 



where Jh is the suitable mesonic operator, projected on the state with 
spatial momentum p. For studying the r\ c and the J/ip, we use the point- 
point operators C75C and 07^ c, respectively. For the vector channel, here we 
consider the transverse component only. Through analytical continuation, 
the Matsubara correlator can be related to the hadronic spectral function 
by an integral equation: 



While the inversion of Eq. (3) to get the spectral function from the 
correlators measured at a finite set of points is clearly an ill-defined prob- 
lem, some progress has been made in the last few years towards tackling 
such problems using Bayesian tools. In the maximum entropy method 8 
the positivity of the spectral function and information about the asymp- 
totic behavior is used to find the most probable spectral function given 
data. Analysis of the hadronic correlators at low temperatures (below T c ) 



&- 1 {p a ,p : T)=pl+p 2 + to 2 + Ren(p, T) + Tmn(p, T) 



(1) 



G H (T,p,T) = (J ff (r,p)4(0, -p)) T 



(2) 
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reveal a complicated structure in the high energy regime of the spectral 
functions, dominated by lattice artifacts 3 . We use this structure as part 
of the asymptotic information in our analysis at higher temperatures 3 . 

Useful information about possible change of state with deconfincment 
can be obtained by comparing the correlators measured above T c with the 
correlators reconstructed from the spectral function below T c . We define 
G rC con,T* (t, p, T) by substituing a(u),p, T) in Eq. (3) with a(w,p, T* where 
T* is usually taken as the smallest temperature below T c available to us. 
This takes into account the trivial temperature dependence of the kernel 
and therefore, any deviation of G(r,p,T) from G rccon .T*(T,p,T) will indi- 
cate a temperature modification of the spectral function 3 . In Fig. 1 the 
results of such a comparison are shown for zero momentum correlators. For 
the pseudoscalar channel, Fig. 1 shows that the reconstructed correlator 
explains the measured correlators completely upto temperatures of 1.5 T c , 
indicating that there is no significant change for these states upto this tem- 
perature. The spectral function extracted from the correlators, shown in 
the right, confirm this. For J/ip, some small deviation of the measured 
correlators is seen from the reconstructed correlators, indicating a possible 
widening 4 or mass shift. Here we added the three (spatial) components of 
the vector, and the behavior is somewhat different from the four compo- 
nent correlator in Ref. [3] due to the change in contribution from electric 
conductivity below and above T c . 




:[fm] «[GeV] 



Figure 1. Comparison of zero momentum correlators at 1.5 T c for vector and pseu- 
doscalar charmonia with those reconstructed from spectral functions at 0.75 T c (left). 
The stars and crosses correspond to a lattice with lattice spacing of 0.02 fm and the 
the open symbols correspond to a lattice with lattice spacing 0.02 fm. (Right) Recon- 
structed spectral functions at 1.5 T c and 0.9 T c for the pseudoscalar channel. Plotted 
is the rescaled spectral function p(uj) = a(ui)/oj 2 . The default model included the peak 
structure at high ui, which is dominated by lattice artifacts. 
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For investigation of the finite momenta charmonia, we use the same 
analysis methods as those at zero momentum described above. Comparison 
of the finite momentum correlators G(r,p,T) at 1.1 and 1.5 T c with the 
reconstructed correlators from a(u,p,T) are presented in Fig. 2. It can be 
seen that with increasing spatial momentum, the correlators above T c show 
larger deviations from the reconstructed correlators, signalling a medium 
effect already at 1.1 T c and increasing with temperature. 
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Figure 2. Comparison of correlators at 1.1 T c and 1.5 T c for pscudoscalar and (trans- 
verse) vector charmonia with those reconstructed from spectral function at 0.75 T c . 



In order to get further insight into the nature of the medium effect, 
we used maximum entropy analysis for a direct reconstruction of a(u>,p, T) 
above T c . As described before, for the reconstruction we use the high energy 
structure of the spectral function obtained at 0.75 T c as part of the prior 
information. The spectral functions at 0.75 T c and 1.1 T c , using such a 
default model, are shown in Fig. 3. Here, in order to avoid any systematics 
due to the small spatial extent and limited number of data points at the 
higher temperature, we use the same physical extent and same number of 
data points both below and above T c , by omitting points at the center for 
the lower temperature. 

Figure 3 shows two interesting features. First, a strong signal from the 
ground state shows up in the spectral function at 1.1 T c even for momenta 
~ 2GeV. Little change in peak strength is seen on crossing deconfmement, 
for any of the momenta. Second, there is an indication of a slower movement 
of the peak position above T c , indicating a modification of the energy- 
momentum dispersion relation. These trends continue at 1.5 T c , with a 
significant bound state showing up with a modified dispersion relation. 

The J /ip moving in the heatbath frame sees more energetic gluons than 
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Figure 3. Spectral functions at 0.75 T c and 1.1 T c for different spatial momenta, for 
lattices with lattice spacings of 0.04 fm. The small horizontal bars indicate the expected 
shift of the peak using relativistic dispersion relation. The same temporal extent has 
been used at both the temperatures. 



the J/ip at rest. The fact that J/ip moving at such momenta show a sta- 
tistically significant peak indicates a stongly nonperturbative nature of the 
plasma. We would like to clarify that we are investigating here J/ip states in 
an equilibriated plasma. This is very different from earlier phenomcnologi- 
cal studies concerned with survival of J / ip with large trasnverse momentum 
from a finite sized plasma, where a large transverse momentum may lead to 
the J I ip spending a smaller time in the plasma, resulting in a larger survival 
probability 9,2 . However, such studies also need to take into account the 
fact that the in-medium dispersion relation for the J/ip will be modified. 
A modified dispersion relation may also explain the large temperature de- 
pendence observed in the screening mass of the J/ip and the r/ c 3 ; however, 
a quantitative statement will require a much more exact determination of 
the dispersion relation than is available at present. 
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